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Abstract. We present new mid-infrared (MIR) diagnos- 
tics to distinguish emission of active galactic nuclei (AGN) 
from that originating in starburst regions. Our method 
uses empirical spectroscopic criteria based on the fact 
that MIR emission from star forming or active galaxies 
arises mostly from HII regions, photo-dissociation regions 
(PDRs) and AGNs. The analysis of the strength of the 
6.2 /im Unidentified Infrared Band (UIB) and the MIR 
continuum shows that UIBs are very faint or absent in 
regions harboring the intense and hard radiation fields of 
AGNs and pure HII regions, where the UIB carriers could 
be destroyed. The MIR signature of AGNs is the presence 
of an important continuum in the 3-10 /im band which 
originates from very hot dust heated by the intense AGN 
radiation field. Using these two distinct spectral proper- 
ties found in our MIR templates, we build diagnostic dia- 
grams which provide quantitative estimates of the AGN, 
PDR and HII region contribution in a given MIR spec- 
trum. This new MIR classification can be used to reveal 
the presence of AGNs highly obscured by large columns 
of dust. 



Key words: Galaxies: active 
frared: galaxies 



Galaxies: starburst - In- 



1. Introduction 

Since the discovery by IRAS of ultraluminous infrared 
galaxies (L(8-1000/xm) > 1O 12 L ) which emit the bulk 
of their energy in the infrared, a large number of stud- 
ies have shown that intense star forming regions as well as 
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AGNs are necessary to explain these high luminosities (see 
ganders fc Mirabel 199£ for a review). Large concentra- 
tions of molecular gas are needed for fueling nuclear star- 
bursts and/or AGNs. Consequently, absorption makes the 
distinction between starburst and AGN activity difficult 
and the estimate of their relative contribution to the total 
infrared luminosity is far from straightforward. Optical 



classification ( 


Veilleux et al. 1995 


Kim et al. 1998 


) is plagued by 



Goldader et al. 1995 



extinction, while ra- 
dio wavelengths, although free from extinction, do not 



result in a definite classification ( 


Condon ct al. 1991 


Lonsdale ct al. 1993, J 


smith ct al. 1998 


i. Based on X-ray 


observations 
ground-based 


(Ogasaka et al. 199 


7, 


Risaliti ct al. 2000), 


infrared observations I 


Roche et al. 1991 


Dudley 199S|, 


Murphy et al. 1999, 


3oifer et al. 200C) and 


on ISO observations (Gcnzel et al. 1998, 


Lutz et al. 1998 


Rlgopoulou et al. 199J 


)), considerable progress has been 



made in defining the fraction of the AGN/Starburst con- 
tributions to the bolometric luminosity of ultraluminous 
infrared galaxies. 

To further examine the AGN and starburst connection, 
we have observed a sample of nearby active and interacting 
galaxies harboring starburst regions and/or AGNs using 
ISOCAM. Due to the unique spectro-imaging capabilities 
of the camera which provide an angular resolution of 4- 
8arcsec between 5 and 16 /jm (200-400 pc at lOMpc), we 
can reveal and study obscured central r egions not vis ible 
at optical wavelengths (A 15fim ~ Ay/70, Mathis 1990 ). 

The goal of this paper is to present an empirical 
method to distinguish and quantify the MIR emission 
coming from starbursts and AGNs. Our approach is based 
on a new MIR diagnostic for estimating the relative im- 
portance of these two main energy sources found in galac- 
tic centers. Our sample and the data reduction methods 
are described in section 2. A brief summary of the typical 
MIR emission encountered in these galaxies is presented 
in section 3, while our MIR templates and diagnostic dia- 
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grams are shown in section 4. We discuss our results and 
their implications in section 5 and the final conclusions 
are given in section 6. 

2. Observations and data reduction 

A large sample of nearby galaxies hosting star formation 
activity and AGN signatu res (see Table |l|) was observed 



with the ISOCAM camera (Cesarsky et al. 1996) on board 
the Infrared Space Observatory (ISO, Kessler et al. 1996f ). 
All observations come from the ISOCAM consortium 
guaranteed time programs and most of them were part of 



e.g. |Charmandaris et al. 1999, 


Mirabel et al. 1998, 


Gallais et al. 1998, 


Mirabel et al. 1999). These galaxies 



have been observed either in the spectro-imaging mode 
with the Circular Variable Filter (CVF) or in the raster 
mode with broad band filters. 

The CVF covers a spectral range from 5 to 16 zim 
with a spectral resolution of ~ 40. For CVF scans, a 
single pointing was made for each galaxy, using 1.5- 
3 arcsec/pixels, giving a total field of view for the 32x32 
pixel array of 0.8-1.6 arcmin. The spatial resolution is 
4-8 arcsec limited by the full width at half maximum 
(FWHM) of the point spread function (PSF). High spatial 
resolution and sensitivity are essential in our analysis in 
order to isolate the central regions in nearby galaxies. Ap- 
proximately 12 exposures of 2-5s each were made at each 
wavelength in addition to ~ 20 exposures at the start of 
each observation to decrease the initial detector transient 
effects. 

The raster maps were made with various broad band 
filters designed to select distinct features such as the 
continuum emission of the very small grains using the 
LW3(12-18 /jm) filter, and the UIBs with the LW2(5- 
8.5 fim) and LW4(5.5-6.5/im) filters. Depending on the 
apparent size of the galaxy, a pixel field of view of 1.5 or 
3 arcsec was selected. Most of the maps were 2x2 rasters 
with 6 pixel overlap consisting of approximately 200 ex- 
posures of 2-5s each. Another 30 exposures were added 
at the start of the observations to decrease the effects of 
transients. 

The ISOCAM data were analyzed with the CAM Inter- 
active Analysis software (CIA0) a nd were calibrated w ith 



the general methods described in [Btarck et al. (1999) . To 



correct for the dark current, a dark model taking into ac- 
count the observing time parameters was subtracted. Cos- 
mic ray hits were removed by applying a multircsolution 
median filtering ( |5tarck et al. 199^ ). Corrections of detec- 
tor memory effects were made with the so-called inversion 
method ( Abergel et al. 1996). The flat field correction was 
performed using the library of calibration data. Finally, in- 
dividual exposures were combined using shift techniques 

1 CIA is a joint development by the ESA astrophysics divi- 
sion and the ISOCAM consortium led by the ISOCAM PI, C. 
Cesarsky, Direction des Sciences de la matiere, C.E.A. France. 



in order to correct the jitter effect due to the satellite mo- 
tions (amplitude ~ 0.5 arcsec). To derive the photometry 
of individual galaxy regions, aperture corrections as well 
as deconvolution techniques ( [Btarck et al. 199*8 ) , were ap- 
plied to account for the overall extension of the PSF in 
a few cases. We estimate that the absolute uncertainty of 
our photometric measurements is ~ 20% while the error in 
the relative uncertainty mainly due to errors on transient 
effect correction is ~10%. 



3. MIR emission of normal, starburst galaxies 
and AGNs. 

MIR spectra (5-16 /im) in galaxies can mainly arise from 
a variety of physical components including: 1) the evolved 
stellar population (Rayleigh- Jeans regime), which can 
dominate in early- type galaxies, 2) emission from the 
ionized interstellar gas, 3) non-thermal emission from 
radio sources and 4) dust particules responsible for the 
underlying continuum and 5) carriers of UIBs centered at 
6.2, 7.7, 8.6, 11.3 and 12.7 ttm. The UIBs are thought to 



Allamandola et al. 1985 



materials ( 


Leger & Puget 1984 


Papoular et al. 1989) and they 



dominate the MIR spectra of gal axies with low or 
intermediate st ar formation activ ity ( Mattila et al. 1999 ) 
as in PDRs (Klein et al. 1999) and diffus e galactic 
regions (|Giard et al. 1988|, |Mattila et al. 19961). The dust 



emitting longward of about 10 fim, is attributed to Very 
Small Grains (here-after VSGs) with radius less than 
10 nm (Desert et al. 199C), and is prominent in regions 
actively forming stars , such as Galactic HII regions (e.g. 
Vcrstractc et al. 1996 , Cesarsky et al. 1996bj ) and in star- 
burst regions. To emit in the MIR, the carriers of the UIBs 
and the VSGs are thought to be stochastically heated by 
the stellar radiation field reaching temperature fluctua- 
tions of the order of 100-1000 K QPuget 
The emission 
through 



Desert et al. 1990) 



Leger 1989 



from ionized gas 
is mainly observed through forbidden lines such 
as [Aril] (6.9 Atm), [NeVI](7.6//m), [ArIII](8.9/Ltm), 
[SiV](10.5imr), [Nell] (12.8 /mi), [NeV](14.3/im) and 
[NeIII](15.6/im). Because of their high ionization poten- 
tial, these lines can be used to trace the hardness of the 
radiation field. The presence of [NeVI] and [NeV] has 
also provided evid en ce for AGNs (|Moorwood et al. 1996 , 
|Gcnzcl ct al. 1998), Sturm et al. 20001) or supernova 



remnants (|Qliva et al. 1999| ). 

In order to understand how the total MIR emission 
in our sample of starbursts and AGNs varies according 
to these different contributions, we first review the well- 
studied MIR properties in resolved nearby galaxies from 
early-type to later type disk galaxies and then address our 
observations in light of these properties. 
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Fig. 1. Typical spectrum of a quiescent star forming 
region. The MIR spectrum (5-16/xm) of a disk re- 
gion of M82 away from the central starburst (thick 
solid line), chosen for its high signal to noise ratio, 
is compared to a modeled spectrum (dashed line) in- 
cluding a continuum and UIBs. The dust continuum 
due to the VSGs is modeled by a black-body curve 
(T=200 K) reddened by a dust absorptio n law with 
Ay =6 mag ( |Dudley fc Wynn- Williams 1997[ ). The UIBs 
marked with the rest wavelength of their peak emission 
are modeled by lorentzian profiles (thin solid lines) (cf. 
Boulanger et al. 1998| ) . Note the faint [Nelll] emission and 
the probable contamination of the [Nell] 12. 8 fim line with 
the 12.7/^m UIB feature. Faint UIBs are also detected at 
5.3, 5.7, 12.0, 13.6, and 14.3 /im as in many other galaxies 
of our sample. 



3.1. MIR emission of "normal" galaxies: early to late 
type 

The MIR spectra of elliptical galaxies, usually poor in cold 
gas and dust, are produced primarily by their evolved stel- 
lar population, and can be shown to resemble a black- 
body continuum with temperatures ranging from 4000 
to 6000K ( [Madden et al. 1997| , [Bosclli et al. 1998| ). This 
stellar emission accounts for most of the MIR emis- 
sion of some elliptical galaxies. However, a non-negligible 
fraction of ellip tical galaxies shows emission from UIBs 
and hot dust (Knapp et al. 1996a , Knapp et al. 1996b , 
Madden et al. 1999 ), as well as emission from non-thermal 
sources. 

The MIR spectra of spiral galaxies, rich in gas and 
dust, are largely dominated by the UIBs. The observed 
UIB features show little spectral variations, and the to- 
tal broad band LW3(12-18/xro)/LW2(5-8.5/mi) flux ra- 



tio is close to 1 (Bosclli et al. 1998, Roussel et al. 1999a, 
Roussel et al. 1999t| ). This ratio has been proposed as 



Fig. 2. MIR spectra showing the variation of the SED 
from a pure HII region, depicted by a solid line (close 
to OB stars in M17, Cesarsky et al. 1996b), to an 
isolated PDR indica ted by a dotted line (NGC 7023, 
|Cesarsky et al. 1996a ). Each spectrum was normalized to 
the continuum at 10 /mi. The evolution is indicated by the 
vertical arrow which shows how the spectrum changes as 
we move from the outskirts of the HII region (dotted line) 
towards its core and near the OB stars (dashed and solid 
lines). Note how the strength of the UIBs progressively 
diminishes and how the continuum emission increases. 



ical" integrated MIR emission of normal spiral galax- 
ies, we present in Figure [j] the spectrum of a region 
in the disk of M82, 45arcsec away from the central 
starburst. The usual UIB features, fitted by lorentzian 
profiles, are prominent. In addition, we notice the 
presence of several faint UIBs at 5.3, 5.7, 12.0, 13.6 



and 14.3/^m (Verstraete et al. 1996, Bcintcma et al. 1996, 
[Sloan et al. 1999 ). Similar MIR spectra are also found 



in regions of the Galactic disk (Cesarsky et al. 1996a, 



Boulanger et al. 1996) 



an indicator of the fraction of the VSG continuum to 
the UIB feature emission. As an example of the "typ- 



3.2. MIR emission of starburst galaxies. 

The MIR spectra of galaxies with evidence of starburst 
activity have distinct features compared to normal spi- 
ral galaxies. Their most prominent characteristic is the 
presence of a very steeply rising continuum at 12-16 /im 
from the VSGs. In starburst regions, this continuum is 
primarily produced in HII regions, tracing regions of mas- 
sive star formation activity. The hard radiation field from 
young stars also excites ionic lines from [Aril], [Arlll], 
[SIV], [Nell] and [Nelll] which can be seen in the MIR 
spectra of starbursts. We display an example of this con- 
tinuum in Figure |^ (solid line) using the spectrum near 
a pure HII region observed in M17, where the MIR emis- 
sion i s almost completely do minated by the VSG contin- 
uum ( |Cesarsky et al. 1996b| ) . The weak intensity of UIBs 
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Table 1. Our galaxy sample 



Source 


RA (J2000) 


DEC (J2000) 


LW2 


LW3 


LW4 


LW3 


LW2 


Spectral type 
















(mJy) 


(mJy) 


(mJy) 


LW2 


LW4 


LW 


CVF 


NGC 1068(Nuc:9")(2) 


02 k 


42 m 


40.6 s 


-00° 00' 


47.8" 


13009 


46616 


10221 


3.58 


1.27 


AGN 


AGN 


Arpll8(Nuc:3")(4) 


02 h 


55 m 


12.2 s 


-00° 11' 


00.8" 


16 


42 


10 


2.62 


1.62 


AGN 


AGN 


NGC3147(Nuc:9")(5) t 


10 h 


16 m 


53.6 s 


73° 24' 


03.3" 


8 


22 


9 


2.60 


0.90 


AGN+ 




Centaurus A(Nuc : 4.5 ")(6) 


13 h 


25 m 


27.6 s 


-43° 01' 


08.8" 


575 


1658 


500 


2.89 


1.15 


AGN 


AGN 


NGC 6814(Nuc:9")(7) t 


19 h 


42 m 


40.6 s 


-10° 19' 


24.6" 


33 


85 


31 


2.59 


1.05 


AGN 






NGC 253(Nuc:7.5")(6) 


00 h 


47 m 


33.1 s 


-25° 


17' 


17.8" 


4703 


15716 


2296 


3.34 


2.05 


HII 


HII 


Arp236(Nuc:4.5")(3) 


0l h 


07 m 


47.5 s 


-17° 


30' 


25.6" 


200 


358 


108 


1.79 


1.84 


PDR 


AGN* 


NGC 1808(Nuc:9")(l) 


05 h 


07 m 


42.3 s 


-37° 


30' 


46.1" 


1074 


1450 


538 


1.35 


2.00 


PDR 


PDR 


M82(Nuc:9")(6) 


09 h 


55 m 


52.2 s 


69° 


40' 


46.9" 


5198 


12720 


2573 


2.45 


2.02 


HII 


HII 


NGC 3256(Nuc:4.5")(l) 


10 h 


27 m 


51.8 s 


-43° 


54' 


08.7" 


442 


1196 


212 


2.70 


2.09 


HII 


HII 


Arp299(A:4.5")(l) 


ll h 


28 m 


34.2 s 


58° 


33' 


46.5" 


325 


1860 


108 


5.73 


3.00 


HII+ 


HII+ 


Arp299(B:4.5")(l) 


n h 


28 m 


31.5 s 


58° 


33' 


40.4" 


505 


1951 


303 


3.86 


1.67 


HII 


HII 


Arp299(C':4.5")(l) 


ll h 


28 m 


31.8 s 


58° 


33' 


49.9" 


76 


232 


36 


3.06 


2.14 


HII 


HII 


Arp299(C:4.5")(l) 


n h 


28 m 


31.2 s 


58° 


33' 


48.9" 


126 


461 


65 


3.66 


1.93 


HII 


HII 


NGC4038(KnotA:6")(8) 


12 h 


01 m 


54.9 s 


-18° 


53' 


12.3" 


23 


135 


13 


5.78 


1.75 


HII 


HII 


Arp 220(9) 


15 h 


34 m 


57.2 s 


23° 


30' 


10.8" 


162 


732 


79 


4.50 


2.07 


HII 


HII 


NGC 6240(9) 


16 h 


52 m 


58.5 s 


02° 


24' 


03.4" 


229 


758 


107 


3.30 


2.14 


HII 


HII 


IRAS 17208-0014(1)* 


l? h 


23 m 


21.9 s 


-00° 


17' 


00.4" 


127 


248 


60 


1.96 


2.11 


PDR 




IRAS 19254-7245(1)+ 


19 h 


31 m 


21.6 s 


-72° 


39' 


20.2" 


111 


264 


84 


2.37 


1.33 


AGN* 




IRAS20551-4250(l) f 


20 h 


58 m 


26.8 s 


-42° 


39' 


00.6" 


123 


425 


62 


3.47 


1.99 


HII 




IRAS 23128-5919(1) 


23 h 


15 m 


47.0 s 


-59° 


03' 


14.0" 


90 


319 


53 


3.54 


1.70 


HII 


HII 




NGC 253(Disk)(l) 


00 k 


47 m 


33.1 s 


-25° 


17' 


17.8" 


339 


681 


171 


2.00 


1.98 


PDR 


PDR 


NGC 520(1) 


0l h 


24 m 


34.8 s 


03° 


47' 


30.8" 


486 


511 


231 


1.05 


2.10 


PDR 


PDR 


NGC 1068(Disk)(2) 


02 h 


42 m 


40.6 s 


-00° 


00' 


47.8" 


246 


296 


137 


1.20 


1.80 


PDR 


PDR 


NGC 1808(Disk)(l) 


05 h 


07 m 


42.3 s 


-37° 


30' 


46.1" 


1615 


2694 


728 


1.67 


2.22 


PDR+ 


PDR 


M82(Disk)(l) 


09 h 


55 m 


52.2 s 


69° 


40' 


46.9" 


1177 


1133 


589 


0.96 


2.00 


PDR 


PDR 


NGC3147(5) t 


10 h 


16 m 


53.6 s 


73° 


24' 


03.3" 


375 


483 


215 


1.29 


1.74 


PDR 




NGC 3256(Disk)(l) 


10 h 


27 m 


51.8 s 


-43° 


54' 


08.7" 


202 


420 


85 


2.07 


2.38 


PDR+ 


PDR+ 


NGC 3263(Nuc:4.5")(l) t 


10 h 


29 m 


13.1 s 


-44° 


07' 


22.0" 


59 


62 


31 


1.03 


1.94 


PDR 




NGC 4676(A)(1) 1 " 


12 h 


46 m 


10.1 s 


30° 


43' 


57.2" 


58 


30 


40 


0.53 


1.44 


PDR+ 




NGC4676(B)(l) t 


12 h 


46 m 


11.4 s 


30° 


43' 


23.1" 


4.11 


2.15 


2.17 


0.52 


1.89 


PDR 




Centaurus A (Disk) (10) 


13 h 


25 m 


27.6 s 


-43° 


01' 


08.8" 


62 


101 


35 


1.62 


1.78 


PDR 


PDR 


NGC6814(6) t 


19' 1 


42 m 


40.6 s 


-10° 


19' 


24.6" 


291 


290 


144 


1.00 


2.02 


PDR 




NGC 7252(1) 


22 h 


20 m 


44.9 s 


-24° 


40' 


41.3" 


142 


185 


67 


1.30 


2.11 


PDR 





Table note: Our galaxy sample is presented in three parts. The top includes nuclei of galaxies containing an AGN, the middle 
contains regions harboring starburst activity and the bottom presents quiescent star forming regions. The spectral type column 
indicates which of the three components (HII, PDR or AGN) provides the dominant contribution. Two values, one according 
to our broad band diagnostic (left column, see also Fig. ^) and one based on the CVF diagnostic (right column, see also Fig.||) 
are presented. The broad band fluxes for most galaxies have been calculated from their CVF spectra. NGC 7252 was observed 
in the CVF mode but due to its weak continuum emission the errors of the CVF diagnostic are large so we include only 
t he LW classification. In NGC 4676(A), the MIR sp ectral energy distribution is probably contaminated by stellar emission 
( [Hibbard fc van Gorkom 1996| , |Read fc Ponman 1998| ) . 

in arcsec, "Disk" = well detected star 
Laurent et al. 1999c|_ 



Notation used: "Nuc"= nucleus with the 
References: 1 This work 



Laurent et al. 2000a; Laurent et al. 1999b 



Laurent et al. 1999. 



diameter 

3 



Charmandaris et al. 2000 



formation 

5 



region m 



irabel et al. 1998 ; Charmandaris et al. 199£ 



the di sk. 
|Laurcnt et al. 1997| ; 6 



Mirabel et al. 1999. The 



identification of specific regions in Arp 299, NGC 4038 and NGC 4676 follows the usual notation found in the literature. Galaxies 
marked with* have been observed only in broad band filter mode (LW). Galaxies marked with + identify those which clearly 
fall outside the limits of our diagnostic, as shown in Figure [] and g and are classified to their nearest MIR spectral type. The 
galaxies marked by * have not been previously classified as AGN-dominant, but have MIR characteristics which our diagnostics 
suggest as those of an AGN. 
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is interpreted as a consequence of the destruction of their 
carriers. However, since the UIBs originate from diffuse 
regions as well as PDRs surrounding HII regions, where 



they peak (Cesarsky et al. 1996a 


3 


Vcrstractc et al. 1996, 


Mattila et al. 1996), 


Iran 1998 


), strong UIBs can be de- 



tected in the MIR spectrum of embedded starburst re- 
gions. Consequently, this depletion of the UIB carriers 
will be further enhanced relative to the observed rel- 
ative increase of the MIR continuum in starburst re- 
gio ns (see spectrum of Knot A in the Antennae Galax- 



Mirabel et al. 1998|) . Therefore, the LW3(12-18 /mi) to 



les, 

LW2(5-8.5 /im) flux ratio increases in strong starburst en- 
vironments. One should note that the overall MIR spectral 
shape can be also considerably affected by the high ex- 
tinction often found in starburst regions. In particular the 
silicate band centered at 9.7/xm can suppress the strength 
of UIBs at 11.3, 8.6 and 6.2 fim relative to that at 7.7/im 



(Lutz et al. 1998, Rigopoulou et al. 1999) as seen in the 



spectrum of Arp220 ( Charmandaris et al. 1999) 



3.3. MIR emission of AGNs. 



Several 
spectra 



studies have 
in AGNs 



already shown that 
also present weak 



MIR 
UIBs 



(Roche et al. 1991, Genzel et al. 1994 Schulz et al. 1998 



Rigopoulou et al. 1999| , |Dudley 1999j ). This effect is 
demonstrated in the spectrum of the nearest AGN of our 
sample located in the radio galaxy Centaurus A (CenA, 
NGC 5128, pVlirabel et al. 1999Q . 

In the MIR observations of CenA, we have suffi- 
cient spatial resolution to disentangle the emission of 
the central regions near the AGN from that due to 
star forming regions of the galactic disk. The absence 
of UIBs in the central 5arcsec (~ 100 pc) possibly indi- 
cates the destruction of their carriers b y the intense UV- 
X-ray radiation f ield from the AGN ( Leger et al. 1989 , 
Allain et al. 1996] ). Moreover, there is a noticeable con- 
tinuum at short wavelengths (3-10 /mi) commonly at- 
tributed to hot dust, associated with the torus of molecu- 



lar gas proposed in the unified model (Pier fc Krolik 1992 



radiation field can heat dust up to evaporation temper- 
atures of ~ 1000K for silicate and ~ 1500K for graphite. 
As a result, the dust continuum emission becomes promi- 
nent at short wavelengths (3-6 /im) in contrast to emission 
from star forming regions which appears at longer wave- 



lengths (e.g. Barvainis 1987). Such a continuum is found 
in the center of all galaxies in our sample known to be 
hosting an AGN such as NGC 1068, NGC 6814 or NGC 



3147 (Laurent et al. 1997). Highly ionized species tracing 
the hard radiation field of the AGN (e.g. [NeV]14.3/tm 
and [NeVI]7.6/im) are detected in NGC 1068. However, 
due to the low spectral resolution of ISOCAM spectra 
(A/AA -40), the [NeV] and [NeVI], if they exist, are 
blen ded with the near by UIBs in all of our galaxies (see 
also |Sturm et al. 2000| ). 

4. MIR diagnostics to distinguish AGNs, 

starburst and quiescent star forming regions. 

Our diagnostics are based on the assumption that the in- 
tegrated MIR emission in galaxies can be represented by 
a sum of contributions originating from (1) regions where 
the dust is predominantly heated by an AGN and (2) re- 
gions where star formation is the main source of energy. As 
mentionned in sections 3.1 and 3.2, star forming regions 
can also be divided in two classes: relatively quiescent re- 
gions where most of the MIR emission originates from the 
PDRs, and starburst-like regions where continuum emis- 
sion at 15 /mi produced in the HII regions is dominant. We 
therefore propose to categorize the spectra from our sam- 
ple using empirical criteria chosen to separate the distinct 
UIB and MIR continuum behaviours. To further exam- 
ine the relative proportion of each contribution, we also 
compare our spectra with three templates, namely, AGNs, 
pure HII regions and isolated PDRs. The relative propor- 
tions of these templates required to reproduce a given MIR 
spectrum allow us to classify the observed region as MIR 
spectral type dominated by AGN, starburst or relatively 
quiescent regions. 



Granato et al. 1997 , Murayama et al. 200Gj) . The alterna- , ^ MIR templates 



tive interpretation of supernova remnants (SNRs) as the 
physical mechanism explaining the AGN phenomenon 
(Terlevich et al. 1992) is not supported by the MIR obser- 
vations. Nearby SNRs, such as Cassiopeia A, Kepler, the 
Crab, RCW 103 and IC443 are all characterized by strong 
ionic and/or molecular line emission and a faint continuum 
(Douvion et al. 1999, Oliva et al. 199£). These character- 
istic features are not found in our integrated MIR spec- 
tra of galactic nuclei which are known to harbor AGNs. 
The MIR spectrum of AGNs is flatter compared to that 
of pure HII regions (see the AGN and HII spectra in 
Figure |). This is in agreement with IRAS observations 
where infrared spectra from active galaxies are generaly 
significantly flatt er with the peak of emission shifted to- 
wards the MIR (|de Grijp et al. 1985[) . Near an AGN the 



To construct the templates of each component, we selected 
three ISOCAM CVF spectra dominated by each of these 
physical processes. These template spectra are presented 
in Figure ||a. 

The strong continuum at 14-15 /im observed in star- 
burst galaxies was modeled by a typical pure HII Galac- 
tic region spectrum from M17 ( |Cesarsky et al. 1996b ). 
UIBs are absent or faint (upper panel of Figure [3|a) and 
the [ArIII](8.9/mi), [SIV](10.5/xm), [Nell] (12.8 /mi) and 
[NeIII](15.6 /xm) are present. A typical starburst spectrum 
can be usually decomposed using the VSG continuum 



with some UIB emission (Fran 1998). Since UIBs origi- 



nate mainly from PDRs or diffuse regions, we used as their 
template, an isolated PDR spectrum from the galactic re- 



6 



Laurent et al.: Mid- infrared diagnostics to distinguish AGNs from starbursts 



flection nebula NGC 7023 flCcsarsky et al. 1996a ) (see the 
lower right part of the Figure |5|a) . 

The AGN component was characterised by the MIR 
spectrum of the central region of CenA (see Fig. || and 
the lower left part of the Figure |a), in which UIBs arc 
absent as in the HII region template. However, in con- 
trast to the HII region template spectrum, the continuum 
at short wavelengths (5-10 /im) is prominent in the AGN 
spectrum. In addition, the equivalent width of emission 
lines of [Nell] and [Nelll] is small, perhaps due to the 
overwhelming continuum contribution. Note that extinc- 
tion affects the AGN template and is greater in CenA 
than in NGC 1068 (see Fig. 

4.2. MIR indicators: [HII] '/[AGN] and [PDR]/[AGN,HII] 

As each MIR template has distinct characteristics, we can 
quantify their relative contribution to the integrated spec- 
trum. We have used two indicators based on the UIB 
strength and the MIR continuum. As an indicator of the 
relative HII region and AGN contribution, we use the ra- 
tio of "warm" (14-15/im) to "hot" (5.1-6.8/xm) continuum 
(Fig. As this ratio decreases, the presence of an AGN 
becomes more evident due to the increase of the relative 
importance of the "hot" continuum. The UIB(6.2/im) to 
continuum (5.1-6.8/im) flux ratio can be used to iden- 
tify the relative contributions of the quiescent star form- 
ing regions to the starburst and/or AGN MIR emission. 
This indicator is very similar to the one proposed by Gen- 
zel which is based on the strength of the 7.7 /jm feature 
( |Gcnzcl et al. 1998 , Lutz et al. 1998). Our estimate of the 
UIB flux has the advantage to be less affected by the 
strong silicate extinction at 9-11 /Ltm while the continuum 
at 5-7 /im is also less contaminated by the VSG contin- 
uum. The effect of a strong absorption on the AGN con- 
tinuum could artificially enhance the presence of a UIB 
at 7.7/im but not at 6.2 /im. Nevertheless, the choice of 
the 6.2 /im which is weaker than the 7.7/im feature ne- 
cessitates spectra with good signal to noise ratio. This 
is the case in our sample of spectra where the UIB at 
6.2 /im is very well dete cted (Figures [§1 and |^ or the spec- 
tra in the Antennae in Mirabel et al. 199q , for some ex- 
amples) . Since our primary goal is to establish the absence 
of UIBs as well as the presence of a well detected AGN 
continuum, at short MIR wavelengths, a criterion based 
on the UIB(5.9-6.8/im)/Continuum(5.T6.8/im) ratio ap- 
plied on high signal to noise spectra gives a more accurate 
diagnostic for revealing the presence of an AGN. As we 
show in the following sections, the relative variation of 
our UIB/Continuum at 6/im is reliable in separating a 
featureless continuum from UIB spectra even if our 6 /im 
continuum measurements may still be slightly affected by 
emission coming from the wing of the 7.7 /tm UIB (see Fig. 



4-3. Diagnostic diagram 

Having defined the three MIR templates, we construct a 
complete library of composite spectra where the contribu- 
tion of each template can vary between and 100%. On 
these composite spectra, we can measure the diagnostic 
ratios presented above and plot them on a 2-dimensional 
diagram (see Fig. |). Three main regions can be identi- 
fied in the corners of this diagram where the signatures of 
AGNs, starbursts, or quiescent star forming regions dom- 
inate the total spectrum. 

The above diagram can be applied only in the galax- 
ies where MIR spectra are available. However, one may 
expand it to cases where only broad band MIR images 
are available. This can be done by using our library of 
spectra to estimate their corresponding fluxes through 
the ISOCAM filters that isolate as closely as possible the 
regions identified in Figure ^. The filters used are the 
LW2(5-8.5/im), LW3(12-18/mi) and LW4(5.5-6.5/mi). As 
we will show in the following sections the broad band di- 
agnostic gives results comparable to those obtained with 
the CVFs. The equivalent fluxes of our CVF templates 
through the three broad band filters are presented in Fig- 
ure [|b. Using the same reasoning as in Figure ^a, we 
estimate the relative strength of UIBs with the LW2(5- 
8.5 /im)/LW4(5. 5-6.5 /im) flux ratio while the ratio of the 
LW3(12-18 /im) /LW2(5-8.5 /im) emission provides an esti- 
mate of the VSG contribution relative to the UIB and/or 
hot dust emission at short wavelengths. Those two diag- 
nostic diagrams have in common the capacity to separate 
individualy regions dominated by AGN, starburst or qui- 
escent star forming regions. The robustness of the two cri- 
teria against observed MIR spectra of galaxies is examined 
in the following section. 

4-4- Application of the diagnostic 



Our MIR indicators were applied to the galaxies of 
our sample (see Table |lj), in which we also included a 
few irregular galaxies such as the blue compact dwarf 
galaxies IIZw40 and NGC 1569 ( |Madden et al. 1999| ) 
and barred spiral ga laxies NGC 1097 and NGC 1365 
(Roussel et al. 1999a) in order to cover the diagram com- 
pletely. In cases where we had adequate spatial resolu- 
tion, we extracted different physical regions in the same 
galaxy. The results are presented in Figure |[ Galaxies 
commonly classified as starbursts, such as Arp220 and 
Arp 299(A), fall toward the top part of the diagram im- 
plying the presence of a large fraction of starburst regions 
compared to quiescent star forming regions. A few nearby 
galaxies with known AGNs are located in the bottom left 
part of the diagram along with our AGN template, the 
nucleus of CenA. Regions of moderate star formation ob- 
served in galactic disks are situated in the region of the 
diagram where we expect to find a dominant PDR contri- 
bution. In Figure 0, we present three typical spectra from 
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Fig. 5. a) Diagnostic diagram based on CVF spectra, b) Diagnostic diagram based on broad band filters. Three distinct 
areas can be defined. In the upper corner, we select spectra dominated by massive HII regions such as those found in 
starbursts. On the lower left, AGN spectra are dominant in a very small region, and finally, PDR spectra fall in the 
lower right part. The solid lines and the associated percentages indicate a constant fraction of one component along 
each line (see section 4.3). 



our sample dominated by UIBs (the disk of M82), an HII 
continuum (the nucleus of M82) as well as an AGN (the 
nucleus of NGC 1068). Based on their position on our di- 
agnostic diagram in Figure q, we estimate the correspond- 
ing fraction of MIR templates according to the Figure [|a. 
This is a first order approximation since we do not take 
explicitly into account the effects of the 9.7 /im silicate 
absorption band on our template spectra. Consequently, 
the observed difference between model and spectra in the 
range of 9 to 11 /im is probably due to the strong silicate 
extinction in embedded starbursts or AGNs. We do not 
attempt to model the emission from the [Nell] (12.8 /im) 
and [NcIII](15.6/im) lines which are independant of the 
dust emission properties. Nevertheless, the overall agree- 
ment presented in Figure ^ is sufficient for modeling the 
general shape of spectra. 

To extend our diagnostics in the case where spectra 
are not available, we present in Figure § the diagnos- 



tic diagram using broad band filters^. As expected, the 
AGN candidates appear in the left part of this diagram. 
For LW3/LW2 ratios close to 1, PDRs dominate the MIR 
emission in our sample as it was the case for normal spi- 
ral galaxies (Bosclli ct al. 1997, Roussel et al. 1999a). For 
LW3/LW2 ratios between 1 and 6, both an AGN and a 
star forming region signature may appear in our spectra 
(e.g. the nucleus of NGC 253 and NGC 1068 in Fig. |). 

We note, however, that the LW2/LW4 flux ratio is well 
adapted to estimate the presence of an AGN contribu- 
tion. For LW2/LW4 ratios lower than 1.5, we clearly see 
that the MIR emission is dominated by the AGN whose 
strong continuum and nearly absent UIBs, contribute to 



2 For galaxies of our sample for which we had only CVF 
spectra, we calculated the equivalent broad band filter fluxes 
taking into account their filter transmission. Note that since 
the CVF ends at 16 /im to match its equivalent LW3 filter to the 
ISOCAM standard LW3(12-18 /im) filter, we have normalized 
the transmission of the latter between 12 and 16 /im. 
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Fig. 6. Diagram based on CVF diagnostics for 33 spectra using the notation of Table 1. The error bars are estimated 
from the rms maps at 1 a for each wavelength along the CVF spectra, added to a systematic error of 10% due to the 
transient correction method. Note that several spiral and blue compact dwarf galaxies included in this figure are not 
listed in the Table [l] since they belong to other samples (Madden ct al. 199E, Rousscl ct al. 1999a). The circles mainly 
located in the bottom left part of the diagram represent the central regions known to contain an AGN. The galaxies 
hosting starburst activity are in the upper part and are marked with diamonds. The spectra of quiescent star forming 
regions are found to lie close to the PDR region and are denoted with squares. The curves demarcate the AGN, PDR 
and HIT spectral type according to Figure |B|a. 



decrease this ratio. For larger LW2/LW4 ratios, emission 
is dominated either by quiescent star forming regions or 
by starburst regions (see Fig. |J). Applying this criterion 
to the nucleus of M82, we clearly see that it is classified 
as a starburst while the nucleus of NGC 1068 has a typi- 
cal AGN signature. Athough this diagram uses only three 
broad band filters, the MIR classification is in good agree- 
ment with results given by the diagnostic diagram based 
on full spectra (see Table [l]). 

Moreover, the diagnostic diagrams of Figures ^| and 
|8| suggest that the well-known ultra-luminous galaxies 
such as IRAS 23128-5919 and NGC 6240 have MIR spec- 
tra which are overwhelmed by a strong starburst signa- 
ture. In the "Super Antennae" (IRAS 19254-7245), which 



harbors a Seyfert 2 nucleus ( Mirabel et al. 199l| ), more 
t han ~70% of its M IR flux originates from the AGN 
( |Laurent et al. 2000b| ). 



One can also use this technique to discover up to 
now unclassified AGN hidden by a large column den- 
sity of dust. Such a candidate is Arp 236 (IC 1623, 
VV114), an infrared luminous system (Lir = 10 1162 
Lq), classified as an early-stage gas-rich merger (M(H2) 
= 5.1xl0 10 M Q , Yun ct al. 1994), and composed of sev- 
eral nuclei and/or compact starburst regions. Previ- 
ous near-infrared and radio studies showed no evidence 
for an AGN ( |Knop et al. 1994 poyon et al. 1995| ). How- 
ever, the MIR sig natures of the eastern region (knot A, 
Doyon et al. 1995 ) place it near the AGN locus in our 
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Fig. 7. Examples of spectra (solid lines) classified as dominated by star forming regions (the disk and the nucleus of 
M82) and by an AGN (the nucleus of NGC 1068) are shown with their corresponding model based on a combination 
of templates (dashed line). For each spectrum, we find the appropriate combination of templates (8% AGN, 50%HII, 
42%PDR for the nucleus of M82; 8%AGN, 15%HII, 77%PDR for the disk of M82 and 77%AGN, 23%HII, 0%PDR 
for the nucleus of NGC 1068) using the location of those sources in the CVF diagram (see Fig. 6). Each composite 
spectrum has been normalised at 14-15 /im. 



CVF diagram (Fig. ^) and, therefore, may contain a non- 
negligible fraction of AGN contribution in its MIR spec- 
trum. Nevertheless, one should note that the possible 
AGN contribution in Arp 236 is not revealed with the LW 
diagnostic (Fig. |J) which is less accurate compared to the 
CVF diagnostic. Knot A is the brightest source in the 
near-infrared. It is also redder and more compact than all 
other sources in the s ystem (not resolv ed at 3.7 fjxa with a 
FWHM of 0.3arcsec, [Knop et al. 1994|) . We estimate that 
~40% of the MIR emission in the region can originate 
from an optically obscured AGN. A MIR spectral classifi- 
cation based on our method is presented in the last column 
of Table [j] for all of our sample. 

5. Discussion 

Our diagnostic diagrams provide a new tool allowing us 
to identify an AGN or star formation (starburst and/or 
quiescent star forming regions) signature in the integrated 
MIR spectrum of a galaxy. 

Although self consistent, our method is based on data 
from a rather diverse sample. To further test its validity, 
we compared it with the MIR diagnostics developed by 
Genzel et al. (1998a) on a complete sample of ULIRGs, 
using ISOPHOT-S. Using the published MIR spectra of 
the brightest sources in the Genzel sample, we calcu- 
lated the strength of the UIB at 6.2 /im in these spectra, 
as well as the equivalent LW2(5-8.5 /im) /LW4(5.5-6.5 /im) 
ratio since this is the principal discriminator between an 
AGN and star forming regions (see Fig. ^|). The classifica- 
tion of galaxies obtained based on our method is in com- 
plete agreement with that of Genzel et al. (1998a). IRAS 
23060+0505, IRAS 19254-7245, Mrk 231 and Mrk 273, 
with an LW2/LW4 ratio of 1.16, 1.36, 1.42 and 1.60 re- 
spectively, are classified by both methods as galaxies con- 
taining a significant AGN contribution in the MIR. Fur- 
thermore, in the wavelength range of ISOCAM (5-16 /im), 
we can better constrain the nature of the continuum at 



7 /im between dust and stellar contribution using the 14- 
15 /im flux, since the LW3(12-18/mr)/LW2(5-8.5/mr) flux 
ratio is greater than 3 for the dust continuum found in 
AGNs or HII regions and less than 0.4 for the stellar con- 
tinuum (Boselli et al. (1997). We are also able to distin- 
guish whether a featureless continuum is due to an AGN 
or a pure HII region (e.g. IIZw40). 

Selective absorption by amorphous silicates, centered 
at 9.7 and 18 /im, can play a crucial role in obscur- 
ing emission over much of the MIR wavelengths. AGNs 
embedded in a large amount of dust could still remain 
undetected by our diagnostic. Studying a sample of 28 
Seyfert 1 and 29 Seyfert 2 with ISOPHOT-S, Schulz et 
al. (1998) have already shown that the high absorption 
in Seyfert 2 galaxies blocks a large fraction (90% on 
the average) of the MIR continuum from the AGN in- 
ner torus. They do detect the silicate feature in emis- 
sion at 9.7 /im which is in favour of a moderately thick 
torus model (Ay^lOO, Granato et al. 1997 ) and rules 
out models with v ery large optical depths (Ay~1000, 
Pier fc Krolik 1992]) . In our sample, among all nuclei (di- 
ameter < 9") which are classified as AGNs based on our 
MIR diagnostics (see Fig. ||), five nuclei are optically se- 
lected as Seyfert 2. Provided sufficient spatial resolution 
is available, we can detect the hot dust continuum asso- 
ciated with the AGN and we are able to disentangle its 
contribution to the MIR emission of the galaxy. If those 
nuclei were observed with the ISOPHOT-S aperture of 
24x24arcsec 2 , we would detect the presence of an AGN 
only in the two nearest Seyfert 2 of our sample NGC 1068, 
CenA and in NGC 6814 which is a Seyfert 1. The dilu- 
tion effect combined with the optical depth of the AGN 
torus would lead to an underestimation of the number 
of detected AGNs. The MIR continuum associated with 
AGNs is not completely suppressed behind the dusty torus 
and may be observed using sufficient spatial resolution 
(~5arcsec) observations on nearby galaxies (D<50Mpc 
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Fig. 8. Diagram based on broad band diagnostics for 45 spectra using the notation of Table 1. The error bars are 
estimated from the rms map at 1 a for CVF spectra or broad band observations added to a systematic error of 10% 
due to the transient correction method. Note that several spiral and blue compact dwarf galaxies not listed in the 
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of quiescent star forming regions lie close to the PDR area (marked with squares). The curves demarcate the AGN, 
PDR and HII spectral type according to Figure ||b. 



with H =75kms~ 1 Mpc" 1 , see |Laurcnt et al. 1999b| ). The 
detection of the AGN continuum in Seyfert 2 is also in bet- 
ter agreement with a dusty torus model producing mod- 
erate absorption (Ay^lOO) as in the model of Granato et 
al. (1997). 

In Figure |l^, we note that the increase of the absorp- 
tion in an AGN spectrum leads to a decrease of the hot 
continuum used in our diagnostic. This further suggests 
that the "true" AGN fraction would be underestimated 
and our detections must be considered as a lower limit. 
An intrinsically low AGN contribution to the MIR spec- 
trum of a source would be difficult to distinguish from the 
star forming emission at 6-7 /im composed by UIBs and 
the VSG continuum. Consequently, we consider that an 



AGN is detected by our diagnostic when the estimated 
fraction of the AGN emission is larger than the contribu- 
tion of starburst or quiescent star forming regions. 

Recently several studies on the origin of the MIR emis- 
sion in galaxy clusters and deep field surveys were per- 
formed with ISOC AM using the LW3 and LW2 broad 
band filters (e.g. Ausscl ct al. 1999, Desert et al. 1999, 



Florcs ct al. 1999] ) . Our diagnostic implies that one cannot 



discriminate between an AGN and a starburst signature 
using just those two filters in the local universe (see Fig. ^|). 
At high redshifts, the two filter bandpasses sample shorter 
rest frame wavelength emission and are less sensitive to 
the VSG continuum. In Figure [fl], we applied a redshift 
correction on our MIR templates at z=0.2 and z=0.4, and 
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Fig. 9. a)In the above figure we display the variation of the AGN-Starburst fraction (indicated in the y-axis) as 
a function of LW3(12-18/im)/LW2(5-8.5 /im) flux ratio based on the three templates of our diagnostic for galaxies 
of the local universe. We divide our library of composite spectra (see Fig. [|) into one component attributed to star 
formation activity (including both the PDR/diffusc and HII region templates) and a second one to AGN. For a given 
LW3 /LW2 ratio and a selected type of star fomation mixture (as indictated by the solid curves of different PDR and 
HII fractions) one can calculate the percentage of the AGN component in the spectrum, b) Same as in a) but for the 
LW2(5-8.5/im)/LW4(5.5-6.5/!m) flux ratio. Note how the use of just the LW3/LW2 ratio is insufficient for detecting 
AGNs. This is not the case for the LW2/LW4 ratio which is relevant to select AGNs from starbursts (see section 4.4). 



we obtained the equivalent diagnostic diagrams derived 
from observations of nearby galaxies. To extrapolate the 
part of the spectrum under 5 /<m, we use a power law for 
the AGN template and a black-body continuum for the 
pure HII region template. We also assume that the PDR 



component under 5 /mi is negligible (Lu et al. 1999). At 
z=0.4, objects with LW3/LW2<3 are classified as AGN- 
dominant in the MIR. For redshifts between 0.5 and 1, 
the LW3 filter samples mainly UIBs (LW2 at z=0), while 
LW2 probes the hot dust emission between 3 and 5 /im 
(assuming that the stellar emission is negligible in spiral 
galaxies and only begins to be substantial below 2/im). 
Consequently, low values of the LW3/LW2 flux ratios in 
these surveys would suggest the presence of an AGN. 

In our diagnostic method, we have implicitly consid- 
ered that the MIR emission from the evolved stellar pop- 
ulation is negligible in spiral galaxies. Even though this 
may be true for galactic disks (Roussel et al. 1999a), ex- 
tra caution is needed for bulges in the center of which 
one finds the AGN. In that case, it is the rising slope of 
the MIR spectrum which would reveal the presence of an 
AGN, even though we cannot rule out completely a pos- 
sible faint contamination from evolved starlight at short 
wavelengths. Moreover, the fact that the core radius of the 
bulge is several orders of magnitude larger than the region 
responsible for the "hot continuum" would in principle fa- 
cilitate their separation in nearby galaxies. For distant 
galaxies where the integrated MIR spectra include the 
whole galactic bulge, it is very difficult to estimate the 



form of the continuum below 6 /im. M83, a well known 
starburst galaxy, is such an example, since we detect con- 
siderable continuum emission at 5-6 /urn (see Fig. ^) . This 
emission is clearly extended outside the unresolved nuclear 
region and can be interpreted as stellar emission from the 
stellar bar (cf. Elmegreen et al. 1998 and Sauvage priv. 
comm. )|. Further more, for distant galaxies the contri- 
bution of the star forming regions surrounding an AGN 
would progressively enter in the beam and dilute the AGN 
MIR signatures. A more detailed discussion on this issue 



is presented in Laurent et al. 1999b 



Our diagnostic method can be further expanded with 
the improved performance of upcoming telescopes, and 
can be used as a guide in scheduling future research pro- 
grams. In addition to an increase in sensitivity, SIRTF will 
provide better wavelength coverage (5-40 /im) than ISO- 
CAM, and better spatial resolution than the ISO-SWS 
observations. In particular, SIRTF would provide a better 
measurement of the continuum for distant faint sources 
and, calculating the depth of both silicate absorption 
bands at 9.7 and 18 /im, a more precise estimate for the ab- 
sorption. An increase in spatial resolution would also per- 
mit a significant improvement in the AGN/Starburst diag- 



3 To further develop our diagnostic for general cases includ- 
ing evolved stellar population, the K band flux at 2.2 /im could 
be used to estimate the contribution of the stellar emission ob- 
served at short MIR wavelengths. A LW2/K or LW3/K flux 
ratio less than 1 would indicate a more important contribution 
from the stellar component. 
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Fig. 11. Variation of the AGN-Starburst fraction at different redshifts (z=0.2 on the left and z=0.4 on the right) as a 
function of LW3(12-18/im)/LW2(5-8.5 /im). We use the same notation as in Figure ^ for the contribution of the star 
formation component. The percentages indicate the fraction of each contribution in the rest-frame wavelengths. The 
LW3 /LW2 ratio which is degenerate for distinguishing AGNs from starbursts at low redshifts can be used for detecting 
the AGN continuum at higher redshifts (i.e. z~0.4) due to the K-correction effect on MIR spectra (see section 5). 



nostics. Both the infrared instruments installed on ground 
based 10 meter telescopes and the Next Generation Space 
Telescope which will cover the MIR spectrum, will probe 
regions of ~10pcs near the AGNs, decreasing significantly 
the effects of beam dilution. 

6. Conclusions 

Using our ISOCAM MIR observations we have obtained a 
new AGN/Starburst diagnostic based on the strength of 
the UIBs at 6.2 /im and the MIR continuum. We conclude 
that: 

1) In AGN spectra (even with a faint starburst con- 
tamination), a strong MIR continuum is present at short 
wavelengths between 3-10 /im. This continuum may be at- 
tributed to very hot dust grains directly heated by the 
central engine. Furthermore, the absence of UIBs suggests 
that their carriers can be destroyed by the strong UV-X 
ray radiation field. 

2) Our MIR diagnostic diagrams can be used to un- 
ravel AGNs that are completely obscured at optical wave- 
lengths. For example, the nucleus of Arp236 which was 
classified as a starburst galaxy triggered by an interaction 
is likely to contain an AGN that contributes ~ 40% to the 
MIR flux. 

3) The emission from ultraluminous infrared galaxies 
can be dominated either by a starburst (e.g. Arp 220) or 
an AGN (e.g. The Super Antennae, IRAS 19254-7245). 
Nevertheless, an AGN can remain partially hidden by the 
torus absorption even in the MIR which may lead to an 
underestimate of its contribution to the whole spectral 
energy distribution. 



4) In the local universe, the LW3(12-18/zm)/LW2(5- 
8.5 /im) ratio alone cannot be used to distinguish AGNs 
from starbursts. However, for galaxies at z=0.4-l, this ra- 
tio is effective in discriminating the AGN from the star- 
burst contribution. 

5) Applications of adapted versions of our diagnostic 
in future instruments such as SIRTF and NGST, which 
will sample the MIR band with higher spatial resolution 
and sensitivity than ISOCAM, should provide better in- 
sights into the heating sources of the interstellar medium 
in galaxies. 
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Fig. 3. Upper le ft panel : A 7 /iin ma p of the dust emission 
in Ccn A (from Mirabel et al. 1999| ). Upper right panel: 
The central region of CenA as seen in the band 7-8.5 /mi 
which traces the most intense UIB emission at 7.7 /tm. Due 
to a good spatial resolution at 7 /im (FWHM = 5 arcsec, 
i.e. 100 pc at 3.5 Mpc), we can separate the nuclear region 
from the disk. Middle panel : The lower curve represents 
the MIR spectrum originating only from the nuclear region 
(100 pc in diameter) whereas the upper curve shows the 
spectrum integrated over a larger region (800 pc in diam- 
eter) which includes both disk structures and the nucleus. 
The silicate absorption cannot be well estimated from the 
global spectrum but is clearly detected in the AGN contin- 
uum. Lower panel : The AGN spectrum (left) represents 
45 % of the energy between 5 and 16 /tm and it contributes 
more than 90 % at wavelengths between 3 and 6 /im. No 
emission associated with star forming regions is detected 
between 3-4 /xm (lower right panel, ISOCAM SW chan- 
nel). While UIBs dominate the MIR emission of the disk, 
they are almost absent in a region of 50 pc radius around 
the AGN. The higher [NelLT] / [Nell] ratio detected in the 
AGN (lower left panel) further indicates the presence of a 
hard radiation field. 
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Fig. 4. The relative strength of the UIB (6. 2 /im) feature 
is estimated by dividing the spectrum between 5.1 and 
6.8 /xm in two different parts. One is the integrated flux 
above the line between 5.9 and 6.8 /im which is attributed 
to the UIB at 6.2 /im and traces essentially the PDRs in 
the star forming regions. The second is the integrated flux 
under the line which is mainly attributed to the AGN con- 
tinuum and/or wing of the UIB at 7.7 /im. The feature to 
continuum indicator is used to estimate the relative con- 
tributions to the heating of the dust from quiescent star 
forming regions compared to the one from an AGN. We 
measure the "hot continuum" (5.1-6.8 /im) which may con- 
tain an AGN contribution and compare it to the "warm 
continuum" (14-15 /im) produced by both star forming re- 
gions and AGNs. As a consequence, the warm to hot con- 
tinuum ratio can be used to quantify the relative contri- 
bution of a starburst continuum and an AGN in the MIR. 
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Fig. 10. The effects of absorption on AGN spectra charac- 
terised by a strong continuum at short wavelength associ- 
ated with an absence of UIBs. In small panels, we display 
the nuclear spectra of NGC 1068 and CenA as well as a 
best fit based on a simple screen model (with the dust ab- 
sorption law of Dudley & Wynn- Williams 1997) applied 



to a power law of spectral index a. The corresponding 
visible absorptions for the fits are 7 and 20 mag respec- 
tively. The large diagram presents how the same power law 
spectrum (a=1.7) affected when absorption ranges from 
Ay=0 to 100. 



